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Abstract 

C^ The crossover between low and high density regimes of exciton-polariton condensates is exam- 

I ined using a BCS wavefunction approach. Our approach is an extension of the BEC-BCS crossover 

Ch theory for excitons, but includes a cavity photon field. The approach can describe both the low 

o 

, O, density limit, where the system can be described as a Bose-Einstein condensate (BEC) of exciton- 



polaritons, and the high density limit, where the system enters a photon dominated regime. In 
contrast to the exciton BEC-BCS crossover where the system approaches an electron-hole plasma, 
the polariton high density limit has strongly correlated electron-hole pairs. At intermediate den- 

cn 

1/-J sities, there is a regime with BCS-like properties, with a peak at non-zero momentum of the 

o 

(^ singlet pair function. We calculate the expected photoluminescence and give several experimental 

L| signatures of the crossover. 

Jh PACS numbers: 71.36.+c,74.78.Na,67.10.-j 




FIG. 1: (a) The occupation probability u| of electron hole pairs of momentum k for a cavity photon 
energy of a; = —1 and photon coupling il = 0.1. Dotted line shows the Is exciton wavefunction 
for comparison. The (b) electron-hole pair, photon, and total particle density, (c) energy per 
electron-hole pair, (d) gap energy versus chemical potential. 



In recent years there is an increasing consensus that a Bose-Einstein condensate (BEC) 
of exciton-polaritons has been realized experimentally [I] [2]. For exciton BECs, it is well 
known that the system crosses over into a BCS state of weakly correlated electrons and holes 
at high density |2l II]- A natural question is then: Does a BEC-BCS crossover also occur 
for exciton-polaritons? Littlewood and co-workers have examined this question and have 
predicted that with increasing density the system transitions from a BEC state to a photon 
BEC state [5]-[8]. In an intermediate density regime under suitable conditions they predict 
a BCS-like regime [6] . The model that they deal with is a model of non-interacting excitons 
coupled to a common photonic cavity. In this model, the excitons do not contain an internal 
electron hole structure and are treated as spins localized on lattice sites. Our purpose here is 
to include the electron and hole components as well as their Coulomb interaction. By doing 
so we find that several new effects are present which have implications on experimentally 
observable quantities. 



The Hainiltonian we consider is 



H = if kin + iicoui + -ffph — fJ-n 
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where e^ and hk are the fermion annihilation operators for electrons and holes with mo- 
mentum k, V{q) = ^ is the Coulomb interaction, nie and nih are the electron and 
hole effective masses, a is the cavity photon annihilation operator, u is the cavity pho- 
ton energy, Q is the coupling strength of the cavity photon to the electrons and holes and 
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is the total particle number operator. We introduce a chemical 
potential /i into the Hamiltonian in order to fix the total particle number. 
We use a BCS wavefunction ansatz of the form 



1$) = cxp[Aat - AV2] Yl U + Vkeihl^ |0) 



(2) 

where M| + f^ = 1 in analogy to the BCS wavefunction used for excitons [Sill]. An additional 
coherent state photonic term is included of the same form as in Ref . [3] . The photon density 
is riph = A^ and the electron-hole density is n^^ = X^fe^fe' giving a total particle density 
n = ?T.pi-j + Ti^^^. 

The BCS wavefunction ^ is equivalent to the solution of the Hamiltonian ([I]) by mean 
field theory [11]. Following the mean field derivation, we assume a pairing Hamiltonian in eq. 
([I]) and restrict fc' = — fc in the attractive Coulomb interaction. In the repulsive Coulomb 
terms we assume a Hartree Fock approximation, and restrict q = k' — k. Hartree terms 
(corresponding to g = 0), are not present due to the assumption of charge neutrality [1]. 
Expanding the operators elck, hj^hk, G\.h_f^, and a around their mean values and keeping 
only linear terms gives 

H = Eo + nD{a + a^) + (to - ^1)0) a + ^ ^fc(e[efc + h\hk) 

k 

- 5^(Afc - n\){eihl^ + h_kek) (3) 



where 

_ ^'k' ^ 

^fc — -j^ ^ — -^fc 

Xfc = ^y(fc-fc')(el,efe.) (4) 

k' 

A, = J2vik-k'){ei,hl,,) (5) 

fc' 

D = J2('W-,) (6) 

Eo = Y, [Xk{eiek) + Akieihl^)] - 2QDX (7) 

fc 

and we have taken aU expectation values as real, 2/m = l/rrie + ^/rnh, A = (a), and 

(e^Cfe) = {h'j^h^k) due to charge neutrality. The photonic part may be diagonalized by 

introducing the operator 5"'' = a''" — c, and demanding that the off-diagonal terms vanish. 

This gives the condition 

A = -i^. (8) 

u — ft 

The remaining part of the Hamiltonian may be diagonalized by a transformation 

Cfc = Mfc7fco + Wfc7fci (9) 

^fc = -^fc7fco + Uklli- (10) 

Demanding that the off-diagonal terms disappear, we obtain the Hamiltonian 

H = eo + iu- fi)B^B + Y,Ek [ll^lk^ + 7li7fci) , (H) 

fc 

where 

eo = XI [^^ ^Ek + Akieihl^) + Xk{eiek) 
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The parameters Uk and Vk satisfy the standard BCS algebra [TT 



UJ — fi 



(ele.)=^^ = ^(l-^) (12) 

(4^-fc) = "fc^fc = — 1^ — (13) 



E, = ^ei + (Afc - fiA)2. (14) 

The gap energy i?^'" is defined to be the value of Ek minimized over all momenta k. In our 
numerical results, we use units such that the momentum is measured in units of l/ag = 



me^/47re/i^, the energy is measured in units Eq = e^/ineaB (SI units used throughout). In 
the low density limit with no photon field, the exciton energy is E/Eq = — 1. Zero detuning 
therefore corresponds to a photon energy of u/Eq = — 1 (the photon energy is negative since 



we measure energies relative to the bandgap energy). The equations (12)-(14) with (l8j) are 
solved self-consistently to obtain our results. 

In Fig. Wbl the occupation probability f | of electron hole pairs is shown. In the low density 
limit (yU ~ —1.4), the distribution coincides with the exciton wavefunction v^ oc 1/(1 + A;^)^ 
[1]. Furthermore, the density of electron- hole pairs and photons is nearly equal in this limit 
(Fig. [lb). Fig. [It shows the energy per electron-hole pair, which approaches an energy of 
E/Eq ^ —1.4. The lowering of the energy is due to the strong coupling and anti-crossing 
of the exciton and photon to form a lower polariton. Finally, the gap energy shown in Fig. 
[lli is equal to the energy of a polariton E'^"' ^ 1.4. The gap energy in this case is the 
energy required to turn a polariton into a free electron hole pair. We thus conclude that 
exciton-polaritons are correctly reproduced in the low density limit. 

As the density is increased. Fig. [T^ reveals that the momentum distribution spreads out 
to higher momentum states. This is precisely the opposite behavior to what is expected in 
the excitonic BEC-BCS crossover. In a standard BCS state, the instability towards forming 
a Cooper pair is weakened with increasing density, because the electron-hole attraction 
becomes increasingly screened by the surrounding electrons and holes. This results in a Vk 
distribution that approaches a Fermi step for high density, for the excitonic case. In the case 
of exciton-polaritons, the electron-hole attraction in fact becomes enhanced with increasing 
density. This is evidenced by the gap energy which increases with density in Fig. [ijl, instead 
of decreasing in the exciton case. 

What is the origin for this enhanced attraction? Fig. [TJd reveals that at high density 
the photon number greatly exceeds that of the electron-hole number. This can be explained 
due to a difference in the particle statistics of the two excitations. Photons are true bosonic 
particles, hence any number of them can be excited with an energy cost uj. Meanwhile, 
electron-hole pairs are fermions, and suffers a phase space filling effect. In order to excite 
more fermionic particles, electrons and holes of increasingly higher momenta and energy 
must be occupied in order to increase the particle number. Thus it is favorable to excite 
photons rather than electron hole pairs to minimize the total energy, explaining the large 
imbalance in these numbers. 




(0 

£ 



2 
1.5 

1 
0.5 



c 


/ 


J 




M=-Eo / 


y 




/ 


X(o=0 



-1.5 



-1 

1^ 



-0.5 




10° 10^ 10^ 10^ 10^ 
n 



FIG. 2: The singlet pair function ^(fc) = uy^V}, for various densities as shown for (a) zero detuning 
Lij = —Eq (b) blue detuning oj = 0. (c) The singlet pair function peak momenta for detunings 
shown, (d) The condensate interaction energy per unit density for various chemical potentials. 

Given that at high density there is inevitably a large number of photons, we may return 
to the original Hamiltonian ([I| to see the consequences. In the high density limit we may 
replace the photon operator by a classical c-number, which we consider to be very large 
a ^ X ^ 1. Discarding all terms that do not contain this factor, we have 



H 



k 



+ AV 



(15) 



The solution of this Hamiltonian is the BCS wavefunction O) with Uk = —Vk = l/"\/2, 
agreeing with the numerical analysis that the average occupation number approaches 1/2 
for all k. The BCS gap in this limit corresponds to a single excitation of the Hamiltonian, 
which has an energy 






2nx. 



(16) 



Thus with increasing density the BCS gap continues to increase in agreement with Fig. [T]l. 

Now we ask to what extent the polariton system possesses BCS features, rather than 

merely a crossover between a polariton BEC state to a photon-dominated regime. The 



criterion given in Ref. [6] is based on a comparison of the energy scales of the BCS gap 
energy i?^'" with the temperature required for condensation to occur ksTj^Ec = ^r/i^n/mpoi, 
where m^^i is the polariton mass. In this definition, the state can be described as "BCS- 
hke" if the energy to disassociate a polariton is lower than the thermal excitation energy 
to prevent a BEC from occurring. In our units, this gives ksTj^Ec/EQ ~ 22000raa^, where 
we used r/ipoi ~ lO^^r/io and as = lOnm, and ttlq is the free electron mass. In terms of Fig. 
flji, this criterion is always much higher than the gap energy -E;^'", for densities exceeding 
n K. 6.5 X 10^ cm^^. Experimentally, such densities have already been achieved, giving the 
result that all current polariton BECs are all in the "BCS-like" regime, according to this 
criterion. 

There is however another sense that the polaritons can be classified as BCS-like. For an 
excitonic BCS state, the singlet pair function '^{k) = u^Vk is peaked near the vicinity of 
the Fermi momentum and has a width of the order of the inverse of the BCS coherence [TO] . 
Figure |2^ shows the singlet pair function for our polariton system, which is peaked at non- 
zero momentum for large densities. The location of this peak has a non-zero momentum 
above a critical density (Fig. ISb). Such behavior is also seen in the excitonic BEC-BCS 
crossover. The difference here is that instead of the singlet pair function becoming sharper 
with increasing density, here the pair wavefunction becomes broader. Using blue-detuned 
(more excitonic) polaritons '^{k) more resembles the excitonic BCS state (Fig. [2b). An 
alternative definition of a polariton BCS phase may be signaled by presence of peak of the 
singlet pair function at non-zero momentum. 

We now turn to the photoluminescence (PL) characteristics of the transition between low 
and high density. Examining the high density limit first, using a similar approximation to 



(15), we use the Hamiltonian 



H = ^Y. Ha + ^fc a^] + oja)a + ^ JZ a^, (17) 

k k 

where a'^ = e^h^f,, tk = h?k'^/m, al = e\.ek = h_,^h^k, and we have exphcitly included 
the semiconductor band gap energy e^ required to create an electron hole pair. The PL 
spectrum is calculated using [12] 

/(e) = -Re / {A^{T)A{0))e-''^^^dT (18) 

TT Jo 

where A is the operator that couples the system to the external PL field. In the case of 
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FIG. 3: (a) The PL intensity of the BCS model in the high density hniit for the photonic PL and 
excitonic PL. A reservoir coupHng 7 is assumed for both cases, giving the Unewidths as shown, (b) 
The photonic PL for the BCS model at all densities. Zero detuning is assumed in all cases. 

polaritons, the PL is generally observed by leakage of the photon field through the micro- 
cavity, hence A = a. It is also possible to observe the PL via a secondary means, from 
the coupling to the exciton field A = a^. This type of coupling is that measured for pure 
excitons and should also be present in principle for polaritons. Experimentally the excitonic 
PL is emitted homogeneously in all directions, whereas the photonic PL is emitted per- 



pendicularly to the sample surface. Evaluating the time correlation function (18) for both 



types of couplings under a mean field approximation, we find a spectrum as shown in Fig. 
[3^. For the exciton coupling, the PL spectrum is identical to the familiar MoUow's triplet 
spectrum found in resonant fluorescence [12] . For the photon coupling only the central peak 



is present. The reason for this difference is illustrated by the single spin version of (17), 



which has eigenstates for high density |±, A^) = {\ ^, N — 1) ±\ ],, N))/\/2, where A^ is the 
number of photons. For large A^, the photon operator does not cause transitions between 
the ± eigenstates: a|±, N) ^ viV|±, A^ — 1). In contrast, the exciton coupling does cause a 
transition cr^|±,A^) ^ (|+,A^— 1) — |— ,A^ — l))/-\/2, giving the side peaks. 

In the low density limit, either type of coupling gives the same PL spectrum, giving the 
familiar emission centered at the lower polariton energy. For the photon coupling case, in 
both the low and high density limits, the action of applying the operator A = a does not 



cause transitions to excited states. Assuming this is true for all intermediate densities, we 
use the results of our ground state BCS wavefunction to obtain the PL emission parameters. 



From the results of the self-consistent equations (12)-(14), we may write down an effective 
Hamiltonian for the BCS theory H^g = ^b'^b + ^Yb'^b'^bb, where 6 is a bosonic operator 
for the effective theory, /^ = ^ and V = \jj^, where A is the sample area and E is an 
energy density. The peak energy of the PL is the energy of adding a single particle to 
the system, which is by definition equal to the chemical potential /i. To determine the 
linewidth of the spectrum, we use the method presented in Ref. [13] to incorporate the 
effect of self-interaction on the PL. The mean field self-interaction energy of the condensate 
nV is shown in Fig. [2ji, which shows that interaction reaches a maximum at intermediate 
density. We attribute this to the fact that at high density the particles are more likely to 
be present in a photonic state, which have no interactions with other particles. Inputting 
the interaction energy into the theory of Ref. [13] gives the PL spectrum Fig. [sb. The 
PL spectrum gradually shifts from the lower polariton energy to the cavity photon energy, 
with a asymmetric linewidth with an exponential tail towards high energy. The linewidth 
decreases again in the high density limit, due to the decreased interactions. 

In the excitonic PL, the side peaks of the Mollow's triplet should disappear at intermediate 
densities, when the saturation effect of the excitons becomes negligible. The central peak 
of the excitonic PL should exhibit a similar behavior to that shown in Fig. [Sjo. We note 
that only the zero center of mass momentum PL is considered in our analysis and we leave 
calculations of dispersion characteristics as future work. 

We have analyzed the crossover between low and high density limits of exciton-polariton 
condensates using a BCS wavefunction approach. Contrary to the exciton case, the electron 
hole pairs have a reduced separation in the high density limit due to the dominant cavity 
photon field. Intuitively we picture this state as a strong cavity photon field continuously 
creating and destroying electron hole pairs at localized positions in the quantum well, result- 
ing in a half occupancy of f^. In the intermediate density regime, the system has BCS-like 
properties in the sense that the pair breaking energy is less than the energy required for de- 
stroying the condensate, and a peak in the singlet pair function develops. The photonic PL 
shifts from the lower polariton energy towards the cavity photon energy with a broadening 
of the linewidth due to the increased interactions. In the high density limit the excitonic 
PL should exhibit a Mollow's triplet type structure, originating from the saturation of the 



electron-hole occupancy. Inclusion of non-equilibrium effects resulting from a coupling to 
an external bath is left as future work. We expect that inclusion of this effect will reduce 
the energy gap due to pair-breaking dephasing processes [H]. Signatures of the high density 
regime should be observable in current experimental systems. 
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